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lineshape  Analysis 
for  Capacitance  Transient  Spectra 


I.  INTRODUCTION 

Capacitance  transient  spectroscopy  has,  in  recent  years,  gained  widespread 
use  in  the  detection  and  study  of  defect  states  in  semiconductor  materials  and 
devices.  1  The  technique  was  originally  developed  by  Lang  and  called  deep-level 

transient  spectroscopy  (DLTS).  Major  variations  on  the  technique  have  been  de- 
3  4 

veloped  by  Kimerling  and  Miller  et  al. 

The  technique  is  based  on  the  work  of  Williams'*  and  Sah  et  al  where  a  tran¬ 
sient  is  produced  in  the  capacitance  of  a  diode  structure,  following  a  sudden  bias 
increase,  and  resulting  from  charge  carrier  emission  from  defects  in  the  depletion 
region.  Determination  of  the  decay-time  constant,  r,  of  the  capacitance  transient 
as  a  function  of  temperature  yields  important  information  about  the  defect  centers. 

(Received  for  publication  1  November  1979) 

1.  Miller,  G.  L. ,  Lang,  D.  V. ,  and  Kimerling,  L.C.  (1977)  Ann.  Rev,  Mat,  Sci. 

1977:377. 

2.  Lang.  D.V.  (1974)  J.  Appl.  Phys.  45:3023. 

3.  Kimerling,  L.C.  (1977)  Radiation  effects  in  semiconductors  -  1976,  Inst,  of 

Physics  Conf.  Series  No,  31,  221. 

4.  Miller,  G.  L. ,  Ramirez,  J.V.,  and  Robinson,  D.A.H.  (1975)  J.  Appl.  Phys. 

46:2638. 

5.  Williams,  Richard  (1966)  J.  Appl.  Phys.  37:341 1. 

6.  Sah,  C.T. ,  Forbes,  L. ,  Rosier,  L.  L. ,  and  Tasch,  Jr.,  A.  F.  (1970) 

Solid-State  Electron.  13:759. 
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In  capacitance  transient  spectroscopy,  periodic  bias  pulses  are  used  to  pro¬ 
duce  a  repetitive  train  of  capacitance  transients.  This  signal  is  then  processed  by 
some  apparatus  which  has  a  peak  response  for  a  particular  t.  As  the  temperature 
of  the  sample  is  varied,  a  peak  in  the  output  of  this  apparatus  will  be  recorded 
each  time  the  emission-time  constant  of  a  defect  matches  that  to  which  the  instru¬ 
ment  is  tuned.  The  lineshape  of  the  spectra  which  are  obtained  depend  on  the 

particular  apparatus  used  to  process  the  capacitance  transient  signal.  Lang  used 

3  4 

a  double  boxcar  integrator.  Kimerling  uses  a  lock-in  amplifier.  Miller*  uses 

4 

an  exponential  correlator.  The  paper  by  Miller  et  al  gives  the  lineshape  for  the 
double  boxcar  and  exponential  correlator,  and  also  for  one  particular  phase  setting 
of  a  lock-in  amplifier. 

Different  phase  settings  are  possible  when  using  a  lock-in  amplifier  due  to  the 
way  it  processes  the  capacitance  transient  signal.  The  tuned  amplifiers  of  the 
input  channel  allow  only  the  fundamental  Fourier  component  of  the  signal  to  pass. 
The  mixer  then  passes  the  signal  unaltered  for  a  half  cycle  and  inverts  the  next 
half  cycle.  It  is  the  phase  of  this  operation  which  is  adjustable.  Finally,  the  out¬ 
put  filters  smooth  this  out  to  a  D.  C.  average  output.  The  phase  of  the  fundamental 
Fourier  component  is  not  constant  in  our  case,  since  as  temperature  is  varied, 
the  time  constant  of  the  exponential  decay  varies,  resulting  in  a  variation  of  the 
phase  of  the  signal  reaching  the  lock-in  mixer. 

If  the  mixer  operation  phase  is  set  to  coincide  with  the  beginning  of  the  capaci¬ 
tance  transient  (end  of  the  bias  pulse)  then  one  gets  the  lineshape  reported  by 
Miller  et  al.^  However,  if  the  phase  is  set  to  be  in  phase  with  the  fundamental 
Fourier  component  at  that  temperature  for  which  it  has  maximum  amplitude,  then 
a  different  lineshape  is  generated.  Some  lock-in  amplifiers  have  a  phase  indepen¬ 
dent  (vector  amplitude)  output  which  gives  still  another  lineshape. 

7  8 

In  previous  work  ’  we  have  used  the  second  mentioned  phase  setting,  choos¬ 
ing  to  be  in  phase  with  the  fundamental  Fourier  component  at  maximum  amplitude 
of  this  component.  This  has  been  shown  to  require  a  phase  setting  of  24.  5°  before 
the  bias  pulse  (beginning  of  transient).  This  phase  setting  is  convenient  since  one 
can  monitor  the  mixer  output  during  a  run  and  verify  that  the  in-phase  point  occurs 
at  the  peak  in  the  output.  This  gives  assurance  that  the  lock-in  amplifier  has  been 
properly  tuned  and  has  not  drifted,  which  can  be  a  concern  at  very  low  frequencies 
(down  to  0.  2  Hz). 

7.  Drevinsky,  P.J.,  Schott,  J.T.,  DeAngelts,  H.M.,  Kirkpatrick,  A.R.,  and 

Minnucci,  J.A.  (1979)  Conf.  Record  of  the  Thirteenth  IEEE  Photovoltaics 
Specialists  Conference  -  1978,  1232.  — 

8.  Schott,  J.T.,  DeAngelis,  H.M.,  and  Drevinsky,  P.J.  (1980)  J.  Electron.  Mat. , 

to  be  published.  - — 

9.  Schott,  J.T.,  DeAngelis,  H.M.,  and  White,  W.R.  (1976)  Air  Force  Technical 

Report,  AFCRL-TR-76-0024.  - ~  "  ~ 
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In  a  previous  report  we  derived  the  expressions  necessary  to  relate  lock-in 
amplifier  frequency  to  emission  time  constant  and  peak  height  to  defect  concen¬ 
tration,  for  this  phase  setting.  Each  time  a  spectrum  at  a  different  frequency  is 
generated,  another  time  constant /temperature  data  point  for  an  Arrhenius  plot 
which  will  yield  the  basic  defect  paramete  ■»  is  received.  However,  this  only  uses 
one  point  of  each  spectrum,  namely,  the  peak  position.  The  rest  of  the  lineshape 
contains  information  about  the  emission  time  constant  at  other  temperatures  as 
well.  In  fact  a  single  spectrum  at  a  single  frequency  contains  all  the  information 
necessary  to  obtain  the  parameters  received  from  an  Arrhenius  plot  using  a  num¬ 
ber  of  spectra  at  different  frequencies. 

In  this  report,  we  derive  the  theoretical  lineshape  for  arbitrary  lock-in 
amplifier  phase,  and  the  expressions  necessary  to  extract  information  from  these 
spectra.  Theoretical  curves  are  compared  with  experimental  spectra.  An  ap¬ 
pendix  gives  a  computer  program  to  generate  theoretical  capacitance  transient 
spectra. 


2.  SIMPLIFIED  LOCK  IN  AMPLIFIER  TREATMENT 

A  simplified  view  of  lock-in  amplifier  operation  is  to  consider  it  to  act  as  a 

correlator  with  a  sine  wave  weighting  function.  This  is  the  view  taken  by  Miller 
4 

et  al  in  comparing  the  lock-in  approach  with  their  exponential  correlator  approach. 
Then  for  a  capacitance  signal  of  the  form  exp  (-t/r),  we  get  a  lock-in  response  of 
the  form 


A 


r/ 


-t/r 

sin 


2  TTt 

L 


dt 


(1) 


2 7r(  1  -  e~L/T) 

L2  .  2 
-Hr  +  4ir 


(2) 


where  the  lock-in  amplifier  is  tuned  to  a  frequency,  F  -  1/L.  If  we  define  a 
dimensionless  parameter,  a,  such  that 


a 


(3) 
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we  then  get  the  lineshape  expression 


Ma ) 


2to2(1  -  e'1/*) 
4i2o2  +  1 


(4) 


Of  course  A  is  implicitly  a  function  of  temperature  T,  through  the  standard  emis¬ 
sion  rate  equation 


a  =  tF 


e^E/kT 

a(v)N 


(5) 


where  k  is  the  Boltzmann  constant,  a  is  the  defect  capture  cross  section,  (v)  is 

the  carrier  thermal  velocity,  and  N  is  the  density  of  states  of  the  appropriate 

4 

band.  This  lineshape  is  the  one  reported  by  Miller  et  al  for  Kimerling's  lock-in 
amplifier  spectra. 


3.  MORE  EXACT  TREATMENT  FOR  0°  PHASE 


For  a  more  careful  treatment,  we  must  consider  the  complete  sequence  of 
input  channel  filtering,  mixing,  and  output  channel  filtering.  The  effect  of  the 
tuned  amplifiers  of  the  input  channel  is  to  pass  only  the  fundamental  Fourier  com- 

q 

ponent  of  the  exponential  train.  In  a  previous  report,  we  showed  this  Fourier 
component  to  be 


Tn 2  7T  t  ,  .  2  71 1 

Fit)  =  aj  cos  +  bj  sin  - 


(6) 


r  ci sin  (Iir+  tan’1T7) 

where 

.  .  2 o(e'1/0,  -  1) 

al - T~5 - 

(4  +  1) 

u  4ira2(e'l/°  -  1) 

1  =  TTZ-r— — 

(4x  o+l) 


(7) 


(8) 


(9) 
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.  4> 


4 


C1 


2g{e~lla  -  1) 

..  2  2  . .  1/2 

(4jr  a  +  1)  ' 


(10) 


(11) 


The  effect  of  the  mixer  operation  and  the  averaging  of  the  output  channel  filters 
can  then  be  given  by 


A 


0 

-  r  I  F(t)  dt 

-L/2 


L'2 

if  F(t)dt 

0 


(12) 


provided  the  phase  of  the  mixer  operation  is  set  to  0  degrees  with  respect  to  the 
beginning  of  the  capacitance  transient.  Substituting  for  the  Fourier  component 
using  Eq.  (6),  we  find  that  the  cosine  term  averages  to  zero,  and  we  are  left  with 


Note  that  this  is  of  the  same  form  as  Eq.  (4)  found  above,  only  differing  by  a  con¬ 
stant  factor. 


4.  LINESHAPE  FOR  ARBITRARY  PHASE 

The  above  lineshape  is  only  correct  for  a  particular  phase  setting  of  the  lock- 
in  amplifier.  As  mentioned  in  the  introduction,  however,  other  phase  settings  are 
sometimes  convenient.  If  we  consider  the  lock-in  mixer  operation  phase  to  be  set 
to  an  arbitrary  angle  ^  (0  <  i)  <  2ir)  before  the  beginning  the  capacitance  transient, 
then  the  lock-in  reponse  will  be  given  by 
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A 


1_ 

L 


-  <t>L,/2ir 

f 


F(t)  dt  -  j- 


(L/2  -0L/2tt) 

/ 


(-L/2 -<t>L/2ir)  -<frL/2ir 


F(t)  dt 


Substituting  for  the  Fourier  component  from  Eq.  (7)  we  have 


(14) 


*  'c‘  fc/ 


-<t>L/2  ir 


.  /  2  7T  t  ^  .  -1  1  \  ,, 

3111  (“IT  +  tan  2i5j  dt 


(-L/2  -4L/2jt) 


(L/2 -<$L/2ir) 

rf 

L/2jt 


in  +  tan'1  jJj)  dt  J 


3i  {‘Tcos  (tan_1  2^  '*)] 


Thus,  substituting  Eq.  (10),  we  have 


4o(l  -  e~^la) 

Ma)  *  y  2  24, 172  c°s 

7r(4  7T  O'  +  1)  ' 


(tan_1  2^’*) 


(15) 


(16) 


We  can  now  consider  some  special  cases.  For  $  =  0,  we  can  construct  a 
right  triangle  with  two  sides  of  length  1  and  2ira,  and  by  calculating  the  hypotenuse, 
immediately  find 

COS  (tan_1  =  (4  AC  l)1/2  U7) 

Making  this  substitution,  we  get  the  lineshape  of  Eq.  (13)  found  above. 

9 

In  a  previous  report  we  found  the  value  of  a,  to  be  denoted  »max,  for  which 
the  amplitude  of  the  fundamental  Fourier  component  (coefficient  c^  given  by 
Eq.  (10)),  is  a  maximum.  To  be  in  phase  with  the  Fourier  component  at  its  max¬ 
imum  amplitude,  we  must  set  the  lock-in  phase  to 


10 


=  tan 


-1 


2  no. 


max 


24.  5° 


The  lineshape  is  then 

Ma)  =  7^2  T~~  ~f^r  coa  f130'1  2^  '  24*5°)  <18) 

jt(4)t  q-  +  1)  '  '  / 

This  is  the  lineshape  appropriate  for  all  our  published  data. 

If  one  has  a  lock-in  amplifier  which  has  a  phase  independent  or  vector  ampli¬ 
tude  output,  then  the  lineshape  is  just  the  general  form  without  the  cosine  factor 


4o(  1  -  e'1/0) 

— 2"S”  .  .1/2 

no  +  1)  ' 


(19) 


Such  instruments  are  really  dual  channel  lock-in  amplifiers  with  one  channel  out¬ 
put  related  to  the  cosine  of  the  phase  angle  as  above,  and  the  other  channel  90°  out 
of  phase  and  related  to  the  sine  of  the  phase  angle.  The  phase  independent  or  vec¬ 
tor  amplitude  output  is  obtained  from  the  sum  of  the  squares  of  these  two  outputs. 

Figure  1  shows  a  plot  of  each  of  the  three  lineshape  equations  plotted  on  a 
temperature  scale  appropriate  for  the  thermocouple  used  in  our  experiments. 

The  curves  are  generated  by  using  Eq.  (5)  to  determine  the  value  of  a  for  each 
temperature  value,  and  then  using  the  appropriate  A(a)  equation  to  determine  the 
amplitude  at  that  temperature.  Figure  2  shows  experimental  spectra  for  two  dif¬ 
ferent  phase  settings.  The  agreement  between  theory  and  experiment  is  excellent. 


5.  DETERMINATION  OF  ACTIVATION  ENERGY  FROM  LINESHAPE 


In  the  introduction,  it  was  pointed  out  that  a  single  spectrum  contains  all  the 
information  necessary  to  determine  the  defect  activation  energy,  provided  more 
than  just  the  peak  maximum  point  is  used.  From  Eq.  (5)  we  can  write 

o  =  C  T-2  eAE/kT  ( 2 C 


-2  1 12 

where  the  T  dependence  comes  from  the  T  '  dependence  of  the  thermal  velocity 
3/2 

and  the  T  '  dependence  of  the  density  of  states  and  where  we  have  assumed  a 
temperature  indpendent  cross  section.  All  other  factors  have  been  lumped  into  the 
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Figure  1.  Theoretical  Lineshapes  for  the  Lock-In  Amplifier  Ver¬ 
sion  of  Capacitance  Transient  Spectroscopy  for  Different  Lock-In 
Amplifier  Phase  Settings  (Mixer  0°  or  24.  5°  Before  Beginning  of 
Capacitance  Transient)  or  for  Phase  Independent  (Vector  Ampli¬ 
tude)  Output.  The  temperature  scale  is  one  appropriate  for 
thermocouple  output 


Figure  2.  Experimental  Spectra  for  the  E-Center  in  Electron  Irrad 
iated  Silicon  for  Two  Lock-In  Amplifier  Phase  Settings 
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constant  C.  We  now  choose  two  easily  definable  points  on  the  lineshape  curve 

A  (a).  For  example,  we  might  choose  a,  and  such  that  A(a.)  =  A(a„)  i  A 

*  *  l  2  in&x 

From  Eq.  (20),  there  are  two  temperatures  corresponding  to  the  two  a's  which 

equally  well  define  these  half-maximum  points.  If  <  a 2  and  <  T2,  then 

corresponds  to  T 2  and  a2  corresponds  to  T^,  since  targe  a  means  small  T.  We 

can  thus  write 

9  AE  'kT, 

“1  -  C  T-*  e  (21) 


and 

,  AE/kT. 

°2  -  C  V  * 

Taking  the  natural  logarithm  of  both  equations,  we  get 
Inuj  fn  C  -  2  In  T2 

ln»2  In  C  -  2  In  Tj 
Subtracting  Eq.  (23)  from  Eq.  (24),  we  get 


AE 


AE 

kT, 


(22) 


(23) 


(24) 


If  we  use  the  approximation 

<nx32 


(25) 


(26) 


(27) 
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for  x  not  too  far  from  unity,  and  If  we  write 


T2  "  Tav  +  * 


T1  “  Tav  *  A 


where  Tav  is  the  average  of  and  Tg,  we 


find 


(A) 

\ti  t2/ 


c*  T  +  =r  9?  T  ST 

av  T  ..  av  max 
av 


Equation  (26)  can  then  be  written 


(28) 


(29) 


AE 


2kT 


max 


(30) 


This  result  has  been  derived  by  others  for  the  0°  phase  lineshape.10  As  mentioned 
above,  one  need  not  use  the  half-maximum  points.  If  one  only  had  half  the  curve, 
for  example,  the  one -half  maximum  point  and  the  full  maximum  point  could  be 
used.  If  the  left  half  of  the  curve  (from  a  temperature  viewpoint)  were  available, 
one  could  get  the  activation  energy  from 


AE 


2kT 


max 


(31) 


Errors  due  to  noise  in  the  experimentally  determined  lineshape  are  just  com¬ 
pounded  when  the  two  points  chosen  are  close  together.  It  is  advisable,  therefore, 
to  use  the  two  half-maximum  points. 

In  order  to  do  all  of  this,  the  appropriate  points  on  the  experimental  curves 
must  be  able  to  be  accurately  determined.  This  means  that  a  good  clean  trace 
with  an  obvious  baseline  must  be  available.  If  there  is  a  weak,  noisy  signal,  or 
one  whose  baseline  is  not  obvious,  it  is  better  to  do  an  Arrhenius  plot  from  a 
number  of  scans. 

10.  Klmerling,  L.  C.,  private  communication. 
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for  the  various  lock-ln  phase  settings. 


Table  1.  Lineshape  Expressions  and  Values  of  Parameters  Defined  in  the  Text, 
for  Various  Lock-In  Amplifier  Phase  Settings 


0° 

24.5° 

Ph-ind. 

A(a) 

4o<l-e-1/“)  ..  -1 

575  ■*) 

4o(l-e‘1/") 

4jt2£»2  +  1 

2  2 ,  ..1/2 cos(tan 
»(4x  a  +1)  ' 

wmxxKUM 

“max 

0.  4243 

0.3485 

0.3485 

A 

max 

0.  1608 

0. 1739 

0.  1739 

sq.  wave 

resp. 

2.  52 

2.33 

2.33 

exp. resp. 

*i 

0.  1291 

0.09657 

0.07559 

"2 

1.9613 

1.  643 

1.775 

f2 

*1 

15.  19 

17.  01 

23.48 

“2  t 

4.622 

4.715 

5.093 

“max 

“max  * 

“1 

3.286 

3.  609 

4.610 

Ratio  of  unit  square  wave  response  (phase  set  for  maximum  response)  to  unit 
exponential  train  response  (phase  set  to  indicated  value). 

^For  calculating  activation  energy  from  left  hand  side  of  peak  lineshape. 

^For  calculating  activation  energy  from  right  hand  side  of  peak  lineshape. 
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Appendix 

Computer  Plotted  Spectra 


The  following  program  plots  capacitance  transient  spectra,  given  the  input  of 
activation  energy,  effective  cross  section  and  amplitude  for  each  level.  The  com¬ 
ment  statements  at  the  beginning  of  the  program  describe  it  more  fully.  It  is 
written  in  Tektest  III  language  for  use  with  a  Tektronix  S-3260  integrated  circuit 
test  system. 

To  run  the  program,  ensure  that  the  proper  disc  (A65B  Hanscom)  is  in  opera¬ 
tion.  Then,  after  typing: 

Ctrl  C 

to  obtain  the  executive  prompter  ($),  type  either: 

DRIVE  0 
or: 

DRIVE  1 

followed  by  a  carriage  return  (CR),  depending  on  which  drive  contains  the  proper 
disc.  Then  type: 


followed  by  a  CR,  and  then: 


IDENT  SCH 


RUN  DLTS 


followed  by  a  CR.  You  then  enter  data  as  requested  by  the  computer. 
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••v-re****^ 


i-mm-  x 


DLTS.EDTSSCH  DISK  NAMES  A65B  HAN8C0M 

DATE  27-APR-79  TIME  10)14  PAGE  1  OF  7 

1.0100  *  TRANSIENT  CAPACITANCE  SPECTROSCOPY 

1.0200  *  WRITTEN  BY 

1.0300  «  JOHN  T.  SCHOTT 

1.0400  *  (LATEST  REVISION-FEB  79) 

1.0500  * 

1.0400  «  THIS  PROGRAM  PLOTS  COMBINED  AND  INDIVIDUAL  TRANSIENT  CAPACITANCE 
1.0700  *  (DLTS)  SPECTRA  AND  SUPERPOSITIONS  OF  THE  TWO.  IT  ASSUMES  THE 
1.0800  •  USE  OF  A  LOCK-IN  AMPLIFIER  WITH  VARIOUS  PHASE  SETTING  OPTIONS. 
1.0900  *  MATERIAL  CONSTANTS  CAN  BE  SET  FOR  EITHER  N-  OR  P-TYPE  SILICON 

1.1000  *  OR  GALLIUM  ARSENIDE.  TEMPERATURE  IS  PLOTTED  ON  A  SCALE  APPRO- 
1.1  100  *  PRIATE  FOR  A  COPPER-CONSTANTAN  THERMOCOUPLE.  HARD  COPY 
1.1200  «  IMAGE  SIZE  MAY  BE  ADJUSTED  TO  DUPLICATE  SCALE  OF  EXPERIMENTAL 
1.1300  *  TRACES  BY  APPROPRIATE  TEMPERATURE  RANGE  INPUT  AND  ADJUSTMENT 
1.1400  *  OF  "SCALE*  VALUES. 

2.1000  SUBROUTINE  DCOORD ( 4 ) > SCOORD ( 4 ) S GRAPHV 

2.2000  SUBROUTINE  VECTRF< 4 > . CURSOR ( 2) . DRAWV( 3) S GRAPHV 
2.3000  ARRAY  E< 10) » SIGMA< 10 ) » AMP < 10) 

2.4000  ARRAY  X<4> . Y<4> »S2<4) 

2.5000  ARRAY  El < 10) .SIGMAI < 10 ) . AMPI < 10 > 

2.4000  ARRAY  EC< 10) .SIGMAC ( 10 ) . AMPC ( 10 ) 

3.0100  PRESET  X— 252 .87. -195.802. 0. 100 
3.0200  PRESET  Y— 4. 197755. -5.535473. Oi 4 .27941 

3.0300  PRESET  S2-.241814E-3. .  1355393E-3. . 9026148E-4.  .720794E-4 

3.0400  CONSTANT  K-8.42E-5 

3.0500  CONSTANT  PHASE- . 4284029 

3.0400  CONSTANT  PI-3.141593 

3.0700  CONSTANT  C1-2E-22 

3.0800  CONSTANT  C2-SE-22 

3.0900  CONSTANT  C3-5E-21 

3.1000  CONSTANT  C4-7E-22 

4.0010  *  ACCEPT  INPUT  INFORMATION 
4.0020  AGAIN  -  0 

4.0030  PRINT  CR 

4.0040  PRINT  "LOCK-IN  AMPLIFIER  PHASE  SETTINGS  ARE  CODED  AS  FOLLOWS) *.CR 

4.0050  PRINT  "  1.  IN  PHASE  WITH  BIAS  PULSE*. CR 

4.0040  PRINT  "  2.  24.5  DEG  BEFORE  BIAS  PULSE". CR 

4.0070  PRINT  "  3.  PHASE  INDEPENDENT  (VECTOR  AMP  OUTPUT) ".CR 

4.0080  PRINT  "ENTER  NUMBER  FOR  PHASE  SETTINGS" 

4.0090  IF( AGAIN)  4.012 

4.0100  ACCEPT  PHAS 

4.0110  GO  TO  4.015 

4.0120  ACCEPT  CHANGE 

4.0130  IF (CHANGE  EG  0)  4.019 

4.0140  PHAS  -  CHANGE 

4.0150  IF (PHAS  NE  INT(PHAS))  4.017 

4.0140  IF( 1<PHAS<3)  4.019 

4.0170  PRINT  "TRY  AGAIN". CR 

4.0180  GO  TO  4.008 

4.0190  PRINT  "MATERIAL  TYPES  ARE  CODED  AS  FOLLOWS) ".CR 
4.0200  PRINT  "  1.  N-TYPE  SILICON*. CR 

4.0210  PRINT  "  2.  P-TYPE  SILICON". CR 

4.0220  PRINT  *  3.  N-TYPE  GALLIUM  ARSENIDE" »CR 

4.0230  PRINT  ■  4.  P-TYPE  GALLIUM  ARSENIDE" »CR 

4.0240  PRINT  "ENTER  NUMBER  FOR  MATERIAL  TYPES" 
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DLTS.EDTSSCH 
DATE  27-APR-79 


time  10:14 


DISK  NAME:  A 65 6  HANSCOM 
PAGE  2  OF  7 


4.0250 
4.0260 
4.0270 
4.0280 
4.0290 
4.0300 
4.0310 
4.0320 
4.0330 
4.0340 
4.0350 
4 . 0360 
4.0370 
4.0380 
4.0390 
4.0400 
4.0410 
4.0420 
4.0430 
4.0440 
4.0450 
4.0460 
4.0470 
4.0480 
4.0490 
4 . 0500 
4.0510 
4.0520 
4 . 0530 
4.0540 
4.0550 
4.0560 
4.0570 
4.0580 
4.0590 
4 . 0600 
4.0610 
4.0620 
4.0630 
4.0640 
4.0650 
4.0660 
4.0670 
4.0680 
4.0690 
4.0700 
4.0710 
4.0720 
4.0730 
4.0740 
4.0750 
4.0760 
4.0770 
4.0780 
4.0790 
4.0800 
4.0810 


IF (AGAIN)  4.028 
ACCEPT  MAT 
GO  TO  4.031 
ACCEPT  CHANGE 
IF (CHANGE  EO  0)  4.045 
MAT  ■  CHANGE 

IF (MAT  NE  INT ( MAT ) )  4.033 

IF ( 1<NAT<4)  4.035 

PRINT  ‘TRY  AGAIN*. CR 

GO  TO  4.024 

IF (MAT  NE  1)  4.03S 

C  •  Cl 

GO  TO  4.045 

IF (MAT  NE  2)  4.041 

C  »  C2 

GO  TO  4.045 

IF (MAT  NE  3)  4.044 

C  «  C3 

GO  TO  4.045 

C  »  C4 

PRINT  *  DO  YOU  WANT  POS.  PEAKS  (ENTER  1).  NEG.  PEAKS  (ENTER 
PRINT  CR.'OR  BOTH  (ENTER  0)?.*  * 

IF ( AGAIN )  4.05 

ACCEPT  YRANGE 

GO  TO  4.056 

ACCEPT  CHANGE 

IF (CHANGE  EQ  0)  4.056 

IF (CHANGE  NE  999)  4.055 

YRANGE-0 

GO  TO  4.056 

YRANGE-CHANGE 

PRINT  'ENTER  TEMP  RANGE  IN  DEG  KELVIN  ( 25<T <350 >  * . CR 

IF( AGAIN )  4.06 

ACCEPT  'TMIN-'.TMIN 

GO  TO  4.063 

ACCEPT  ' THIN- * .CHANGE 

IF (CHANGE  EO  0)  4.069 

TM IN-CHANGE 

IF ( 25<TMIN<350 )  4.066 

PRINT  'TEMP  OUT  OF  RANGE  25<T<350*.CR 

GO  TO  4.057 

IF ( AGAIN )  4.069 

ACCEPT  *  TMAX- ' . TMAX 

GO  TO  4.072 

ACCEPT  * TMAX-* .CHANGE 

IF ( CHANGE  EO  0)  4.081 

TMAX-CHANGE 

IF (TMAX  LE  350)4.075 

PRINT  'TEMP  OUT  OF  RANGE  25<T<350*.CR 

GO  TO  4.066 

IF ( TMAX  GT  THIN)  4.078 

PRINT  'TRY  AGAIN*. CR 

GO  TO  4.056 

IF ( AGAIN)  4.081 

ACCEPT  'WHAT  IS  THE  LOCK-IN  FREQUENCY?:  *,F 
GO  TO  4.086 

ACCEF'T  *«HAT  IS  THE  LOCK-IN  FREQUENCY?:  ’.CHANGE 
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OLTS • EDT ! SCH 
HATE  27- APR-79 


TIME  10114 


DISK  NAME!  A65B  HANSCOM 
c-AGE  3  OF  7 


4,0820 
4.0830 
4.0840 
4.0850 
4.0860 
4.0870 
4.0880 
4.0890 
4.0900 
4.0910 
4.0920 
4.0930 
4.0940 
4 .0950 
4.0960 
4.0970 
4.0980 
4.0990 
4.1000 
4.1010 
4.1020 
4.1030 
4.1040 
4.1050 
4.1060 
4.1070 
4.1080 
4.1090 
4.1100 
4.1110 
4.1120 
4.1130 
4.1140 
4.1150 
4.1160 
4.1170 
4.1180 
4.1190 
4.1200 
4.1210 
4.1220 
4,1230 
4.1240 
4.1250 
4.1260 
4.1270 
4 . 1 280 
4,1290 
4.1300 
4.1310 
4.1 320 
4 . 1 330 
4 . 1 340 
4.1350 
4.1360 
4.1370 
4  .  t  3!|(, 


IF (CHANGE  EG  0)  4.088 
F “CHANGE 
GO  TO  4.088 
IP (AGAIN)  4.088 

ACCEPT  ' HOW  MANY  PEAKS  FOR  THE  COMBINED  SPECTRUM?!  ',N 
GO  TO  4.094 

ACCEPT  ' HOW  MANY  PEAKS  FOR  THE  COMBINED  SPECTRUM?!  *  r CHANGE 
IF ( CHANGE  £Q  0)  4.094 
IF ( CHANGE  NE  999)  4.093 
N=0 

GO  TO  4.124 
N=CHANGE 

I F ( 1 <N< 1 0  >  4.124,4.097,4.095 
PRINT  'THAT'S  A  BIT  MUCH' ,CR 
GO  TO  4.085 
LOOP  4.123  J=1 , N 
PRINT  CR 

PRINT  'ENTER  ACTIVATION  ENERGY  IN  EV  FOR  PEAK  *,J!I0C,'!  ' 

IF ( AGAIN )  4.103 
ACCEPT  EC ( J ) 

GO  TO  4.106 
ACCEPT  CHANGE 
IF ( CHANGE  EG  0)  4.106 
EC  <  J ) =CHANGE 

PRINT  'ENTER  CROSS  SECTION  IN  SG  CM  FOP  PEAK  ',J!I0C,‘!  ' 

IF ( AGAIN )  4.11 
ACCEPT  SIGMAC(J) 

GO  TO  4,113 
ACCEPT  CHANGE 
IF ( CHANGE  EG  0)  4.113 
SIGMAC(J) =CH ANGE 

PRINT  'ENTER  AMPLITUDE  (POS.  OR  NEC.)  IN  CM  FOR  PEAK  '.JIIOC,'.'  ' 
IF (AGAIN)  4.117 
ACCEPT  AMPC(J) 

GO  TO  4.123 
ACCEPT  CHANGE 
IF (CHANGE  EG  0)  4.123 
IF(CHANGE  NE  999)  4.122 
AMF’C  (  J )  =0 
GO  ro  4.123 
AMFC ( J ) “CHANGE 
CONTINUE 

PRINT  CR  >  *  DO  YOU  WANT  ANY  PEAKS  PLOTTED  INDIVIDUALLY?* ,CR 
IF  (AGAIN)  4.128 

ACCEPT  'IF  SO,  ENTER  THE  NUMBER!  OTHERWISE  HIT  CR!  *,NI 
GO  TO  4. 134 

ACCEPT  'IF  SO,  ENTER  THE  NUMBER!  OTHERWISE  ENTER  999!  ' , CHANGE 
IF (CHANGE  EG  0)  4.134 
IF ( CHANGE  NE  999)  4.133 
NI=0 

GO  TO  5.01 
N I =CHANGE 

IF1 1  N I  10)  ^.Ol ,4 . 137,4. 135 

PRINT  'THAT'S  A  BIT  MUCH.'.CR 

GO  TO  4,125 

IF (N  NE  NI '  4.154 

PRINT  CR 
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ut  ut  u»  ui  in  ui  m  w  ci  m  ui «  cn 


dlts.edtisch 

DATE  27-APR-79 


TIME  10S14 


DISK  NAME:  A65B  HANSCClM 
PAGE  4  OF  7 


4.1390 
4.1400 
4.1410 
4. 1420 
4.1430 
4.1440 
4.1450 
4.1460 
4.1470 
4.1480 
4.1490 
4.1500 
4.1510 
4.1520 
4.1530 
4.1540 
4.1550 
4.1560 
4.1570 
4.1580 
4.1590 
4.1600 
4.1610 
4.1620 
4.1630 
4.1640 
4.1650 
4.1660 
4.1670 
4.1680 
4.1690 
4.1700 
4.1710 
4.1720 
4.1730 
4.1740 
4.1750 
4.1760 
4.1770 
4.1780 
4.1790 
4.1800 
4.1810 

.0100 
.0200 
.0300 
.0400 
.0500 
.0600 
.0700 
.0800 
.0900 
.  1000 
.1100 
.1200 
.1300 


PRINT  ’IF  THESE  ARE  THE  SAME  PEAKS  AS  FOR  THE  COMBINED  SPECTRUM 
IF(AGAIN)  4.143 

ACCEPT  'ENTER  +li  OTHERWISE »  HIT  CR:  '.SAME 
GO  TO  4.147 

ACCEPT  'ENTER  +1?  OTHERWISE.  ENTER  999!  '.CHANGE 

IF (CHANGE  EQ  0)  4.147 

IF (CHANGE  EQ  999)  4.155 

SAME=CHANGE 

IF ( SAME  EQ  0)  4.155 

LOOP  4.152  J-l.N 

El ( J) =EC( J) 

SIGMA I ( J  > “SIGMAC ( J ) 

AMPI(J)=AMPC(J> 

CONTINUE 
GO  TO  5.01 
SAME=0 

LOOP  4.181  J=1.NI 
PRINT  CR 

PRINT  'ENTER  ACTIVATION  ENERGY  IN  EV  FOR  PEAK  I-'.JSIOC.':  * 

IF ( AGAIN )  4.161 
ACCEPT  EI(J) 

GO  TO  4.164 
ACCEPT  CHANGE 
IF (CHANGE  EQ  0)  4.164 
El ( J ) =CHANGE 

PRINT  'ENTER  CROSS  SECTION  IN  SQ  CM  FOR  PEAK  I-'.JSIOC. '8  * 

IF( AGAIN)  4.168 
ACCEPT  SIGMAI(J) 

GO  TO  4.171 
ACCEPT  CHANGE 
IF ( CHANGE  EQ  0)  4.171 
SIGMAI ( J ) “CHANGE 

PRINT  'ENTER  AMPLITUDE  (POS.  OR  NEG . )  IN  CM  FOR  PEAK  I-'.JIIOC. 
IF(AGAIN)  4.175 
ACCEPT  AMP I ( J ) 

GO  TO  4.181 

ACCEPT  CHANGE 

IF  <  CHANGE  EQ  0)  4.161 

IF  (CHANGE  NE  999)  4.18 

AMF'I  (  J )  3=0 

GO  TO  4.181 

AMF'I  (  J )  ^CHANGE 

CONTINUE 

DRAW  AXES 

PRINT  ERASE 

VECTRF (40.760.40.110) 

VECTRF ( 7 . 700  >12.716) 

VECTRF ( 12.716.17.700) 

VECTRF (17, 700. 7. 700) 

CURSOR (21 ,700) 

PRINT  'C* 

VECTRF (40. 110, 1000, 110) 

VECTRF ( 1000.110,1000,780) 

VECTRF ( 1000,780,40.780) 

CURSOR (967, 66) 

PRINT  TEMP* 


DLTS • EDT ! SCH 
DATE  27-APR-79 


TIME  10!14 


DISK  NAME:  A65B  HANSCOM 
PAGE  5  OF  7 


6.0100 
6.0200 
6.0300 
6.0400 
6 . 0500 
6.0600 
6.0700 
6.0800 
6.0900 
6.1000 
6.1100 
6.1200 

6.1300 

6.1400 

6.1500 
6.1600 

6.1700 
6.1800 

6.1900 
6.2000 
6.2100 
6.2200 

7.0100 

7.0200 

7.0300 

7.0400 

7.0500 

7.0600 

7.0700 

7.0800 

7.0900 

7.1000 

7.1100 

7.1200 

7.1300 

7.1400 

7.1500 
7.1600 

7.1700 
7.1800 

7.1900 
7.2000 
7.2100 
7.2200 
7.2300 
7.2400 
7.2500 
7.2600 
7.2700 
7.2800 
7.2900 
7.3000 
7.3100 
7.3200 
7.3300 
7.3400 


DRAW  SCALE  ON  HORIZONTAL  AXIS 
TMIN=INT ( TMIN  > 

TMAX=INT (TMAX) 

TK=TMAX 
CALL  9.01 
EMFMAX=EMF 
TK=TMIN 
CALL  9.01 
EMFMIN=EMF 

SCAL£=(EMFMAX-EMFMIN>/70 

DCOORD  <  EMFMI N- 1 . 54SCALE . EMFMAX , -22 , 4 ) 

SCOORDC 19, 1000.88. 114) 

LOOP  6.21  TK=TMINiTMAX 
IF< TK/50  NE  INT (TK/50) )6.2 
CALL  9.01 
DRAUV  C  EMF  .4,0) 

DRAWV  C  EMF  * -4 1 1 ) 

DRAWV <£MF-1.S*SCALE, -22,-1 > 

PRINT  tk:ioc 

CONTINUE 

CONTINUE 

DRAWV (EMFMIN .0,0) 

CALCULATE  AND  DRAW  DLTS  SPECTRA 

(‘SCALE*  VALUES  MUST  BE  ADJUSTED  FOR  HARD  COPY  IMAGE  SIZE) 

IF ( PHAS  NE  1)  7.06 

SCALE  =  2.31 

GO  TO  7.07 

SCALE  =2.48 

SCOORD< 40. 1000. 110.780) 

IF  <  YRANGE )  7.09.7.11,7.14 
DCOORD  <  EMFMIN, EMFMAX , -SCALE » 0 ) 

GO  TO  7.15 

DCOORD  <  EMFMIN , EMFMAX . -SCALE/2 , SCALE/2 ) 

VECTRF ( 40 , 445 . 1000.445 ) 

GO  TO  7.15 

DCOORD ( EMFMIN . EMFMAX »  0 , SCALE ) 

IF (N  LT  1)  7.31 
-COMBINED  SPECTRUM 
LOOP  7.21  J=1»N 
E< J)=EC< J) 

SIGMA ( J ) =SIGMAC ( J  > 

AMP< J)=AMPC( J) 

CONTINUE 

LOOP  7.3  TK=TMIN» TMAX » . 5 
ASUM=0 

LOOP  7.27  J=1,N 

CALL  10.01 

ASUM  =  ASIIM  +  A 

CONTINUE 

CALL  9.01 

DRAWV ( EMF  »  ASUM » 1 ) 

CONTINUE 
IF <NI  LT  1)  8.01 
-INDIVIDUAL  SPECTRA 
LOOP  7.37  J=1,NI 
E  <  J ) =EI <  J ) 
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DLTS.EDTJSCH  DISK  NAME •  A65B  HANSCOM 

DATE  27-APR-79  TIME  10!  14  PAGE  6  OF  7 

7.3500  SIGMA< J)=SIGMAI< J) 

7.3600  AMP  <  J ) = AMP I  <  J ) 

7.3700  CONTINUE 

7.3800  LOOP  7.49  J=1.NI 

7.3900  DRAWV<EMFMIN.O.O) 

7.4000  LOOP  7.48  TK=TMIN. TMAX# . 5 
7.4100  CALL  10.01 

7.4200  CALL  9.01 

7.4300  IF  <  TK  EQ  INT(TK))  7.46 

7.4400  DRAWV<EMF.A. 1 > 

7.4500  GO  TO  7.47 

7.4600  DRAUV ( EMF  >  A . 0 ) 

7.4700  CONTINUE 

7.4800  CONTINUE 

7.4900  CONTINUE 

8.0100  *  PRINT  INPUT  INFO 
8.0200  CURSOR ( 770 »  44 ) 

8.0300  IF (MAT  NE  1)  8.06 

8.0400  PRINT  *N-SI* 

8.0500  GO  TO  8.13 

8.0600  IF  < MAT  NE  2)  8.09 

8.0700  PRINT  *P-SI* 

8.0800  GO  TO  8.13 

8.0900  IF<MAT  NE  3)  8.12 

8.1000  PRINT  •N-GAAS* 

8.1100  GO  TO  8.13 

8.1200  PRINT  'P-GAAS' 

8.1300  IF < PH AS  NE  1)  8.16 

8.1400  PRINT  •  (0  DEG)- 

8.1500  GO  TO  8.2 

8.1600  IF < PH AS  NE  2)  8.19 

8.1700  PRINT  *  <24.5  DEG)- 

8.1800  GO  TO  8.2 

8.1900  PRINT  -  (PH  IND >  * 

8.2000  CURSOR <714 i 22) 

8.2100  PRINT  FJRl.'HZ- 

8.2200  CURSOR <742.0) 

8.2300  PRINT  TMIN:iO.‘<  T  <  ’.TMAXJIOC 
8.2400  CURS0R<0 »  44 ) 

8.2500  IF ( N  EQ  0)  8.3 

8.2600  LOOP  8.29  J=1,N 

8.2700  PRINT  ’PEAK  - . J ! IOC . - ! ~H " . EC ( J ) ! R3C . -  EU!  - . SIGMAC < J > 1 E2C 

8.2800  PRINT  *  SQ  CM!  AMPL= - . AMPC ( J ) J F3C . -CM - » CR 

8.2900  CONTINUE 

8.3000  IF  <  NI  EQ  0)  8.36 

8.3100  IF < SAME  EQ  1)  8.36 

8.3200  LOOP  8.35  J=1»NI 

8.3300  PRINT  'PEAK  I - - , J! IOC. * t "H - . El < J ) ! R3C r "  EU)  • » SIGMAI <  J >  J E2C 
8.3400  PRINT  -  SQ  C Mi  AMPL=  *  >  AMF’I  ( J) I F3C .  *CM - » CR 

8.3500  CONTINUE 

8.3600  PRINT  CR.CR 

8.3700  GO  TO  11.01 

9.0100  *  CALCULATION  OF  APPROXIMATE!'  EMF  AS  FUNCTION  OF  T 

9.0200  *  (SEE  'TEMPERATURE.  ITS  MEAS  $  CTRL  IN  SCI  1  IND.-  UOL  4,  PT  3 

9.0300  *  ED.  H.H. PLUMB  ( PITTSBURG . 1972 ) ,  PG.  1613.) 
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9.0400  TC=TK-273 

9.0500  1=2 

9.0600  IF(TC  LT  X<3))9.09 

9.0700  1=3 

9.0800  GO  TO  9.13 

9.0900  IF  <  TC-X< I ) )  9.12*9.13*9.1 

9.1000  1=1+1 

9.1100  GO  TO  9.09 

9.1200  1=1-1 

9.1300  HT1«TC-X(I) 

9.1400  HT2=TC-X< 1+1 > 

9.1500  PR0D=HT1*HT2 
9.1600  H=X ( 1+1 ) -X ( I ) 

9.1700  DELY=<  Y < 1+1 > -Y< I >  >/H 

9.1800  S3=(S2( 1+1 )-S2( I ) )/H 

9.1900  SS2=S2< I )+HTl*S3 

9.2000  DELSQS=(S2( I >+S2(I+l >+SS2>/6 
9.2100  EMF=Y(I)+HT1*DELY+PROD*DELSOS 
9.2200  RETURN 

10.0100  *  CALCULATION  OF  DLTS  CURVE  AMPLITUDE 

10.0200  *  (  *  IF '  STATEMENTS  AVOID  EXCEEDING  RANGE  OF  10**-38  TO  10**+38> 

10.0300  IF(E( J)/(K*TK>  LT  86)  10.06 

10.0400  A=0 

10.0500  GO  TO  10.17 

10.0600  P0WR1=EXP(E( J)/(K*TK) ) 

1 0 . 0700  ALPHA=F*C*F'0WR1 / (SIGMA ( J  >  *TK**2  > 

10.0800  IF (ALPHA  GT  .013)  10.11 

10.0900  F0WR2-0 

10.1000  GO  TO  10.15 

10.1100  P0UR2=EXP( -1/ALPHA) 

10.1200  IF(E(J)/(K»TK)  LT  43)  10.15 
10.1300  A1=AMP< J)*2#(1-P0WR2)/PI**2 
10.1400  GO  TO  10.19 

10.1500  IF( PH AS  NE  1)  10.18 

10.1600  A  =  AMP( J)*8*ALPHA**2*( 1-P0UR2)/(4*PI**2*ALPHA**2+1 > 

10.1700  GO  TO  10.23 

10.1800  A1  =  AMP(J**4*ALPHA*(l-P0WR2)/(PI*SQRT(4*PI»*2*At,FHA**2+l>> 

10.1900  IF ( PHAS  NE  2)  10.22 

10.2000  A  =  A1*C0S(ATAN(1/(2*PI*ALPHA> >-PHASE> 

10.2100  GO  TO  10.23 
10.2200  A  =  A1 
10.2300  RETURN 

11.0100  *  ALLOW  REPEAT  WITH  SOME  INPUT  INFO  CHANGED 

11.0200  PRINT  "IF  YOU  WANT  A  REPEAT  WITH  SOME  INPUT  DATA  CHANGED  *  '  *  CR 
11.0300  ACCEPT  'HIT  CR*  OTHERWISE  ENTER  IS'.  REPT 
11.0400  IF (REPT  NE  0)  11.11 

11.0500  PRINT  CR  *  *  AS  INPUT  DATA  IS  REQUESTED*  IF  YOU  WISH  THAT  QUANTITY*  *CR 
11.0600  PRINT  'UNCHANGED*  SIMPLY  HIT  CR.'.CR 

11.0700  PRINT  '(NOTES  TO  CHANGE  A  NON-ZERO  VALUE  TO  ZERO.  YOU  MUST'.CR 
11.0800  PRINT  'ENTER  999.)'.CR 
11.0900  AGAIN-REPT+1 

11.1000  GO  TO  4.003 

11.1100  STOP 
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